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Abstract. The control of magnetic properties by means of an electric field 
is an important aspect in magnetism and magnetoelectronics. We here utilize 
magnetoelastic coupling in ferromagnetic/piezoelectric hybrids to realize a voltage 
control of magnetization orientation at room temperature. The samples consist of 
polycrystalline nickel thin films evaporated onto piezoelectric actuators. The magnetic 
properties of these multifunctional hybrids are investigated at room temperature as 
a function of the voltage controlled stress exerted by the actuator on the Ni film. 
Ferromagnetic resonance spectroscopy shows that the magnetic easy axis in the Ni 
film plane is rotated by 90° upon changing the polarity of the voltage V p applied to 
the actuator. In other words, the in-plane uniaxial magnetic anisotropy of the Ni 
film can be inverted via the application of an appropriate voltage V p . Using SQUID 
magnetometry, the evolution of the magnetization vector is recorded as a function of 
V p and of the external magnetic field. Changing V p allows to reversibly adjust the 
magnetization orientation in the Ni film plane within a range of approximately 70°. 
All magnetometry data can be quantitatively understood in terms of the magnetic 
free energy determined from the ferromagnetic resonance experiments. These results 
demonstrate that magnetoelastic coupling in hybrid structures indeed is a viable option 
to control magnetization orientation in technologically relevant ferromagnetic thin films 
at room temperature. 
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1. Introduction 

Multifunctional material systems unite different, e.g., electric and magnetic, 
functionalities in either a single phase or a heterostructure. They thus are of great 
fundamental and technological interest. Given that the electric and magnetic properties 
furthermore are coupled, it becomes possible to control either the magnetization by 
the application of an electric field or the electric polarization via a magnetic field 
alone. An electric field control of magnetization is particularly appealing, as it removes 
the need to generate magnetic fields of sufficient strength for magnetization switching 
on small length scales - thus enabling novel concepts for high density magnetic data 
storage applications. Thus, techniques to control the magnetization by means of electric 
fields or currents have recently been vigorously investigated and several schemes for 
an electric control of magnetization have been reported. These include the spin-torque 
effect [U El El H] in spin- valves and magnetic tunnel junctions and the direct electric field 
control of magnetization in intrinsically multiferroic materials [S, El d El El QUI CCD, CE2] 
or in ferromagnetic/ferroelectric heterostructures [TBI [14]. A third, very attractive 
approach for the electric field control of magnetization takes advantage of the elastic 
channel, i.e., magnetostrictive coupling [IS], [16l [17l HEl UHl EDI EH [22], [23j . We here show 
that magnetoelasticity - that is the effect of lattice strain on magnetic anisotropy - 
in a polycrystalline nickel film/piezoelectric actuator hybrid structure allows to switch 
the magnetic easy axis in the Ni film by 90° at room temperature by simply changing 
the polarity of the voltage applied to the actuator. This can be used to achieve either 
an irreversible magnetization orientation control of 180° or a reversible control close 
to 90°, depending on whether the magnetization orientation is prepared in a local or 
global minimum of the free energy via a magnetic field sweep prior to the electric 
control experiment. We use ferromagnetic resonance (FMR) to quantify the effect of an 
electric field on magnetic anisotropy and superconducting quantum interference device 
(SQUID) magnetometry to directly record the evolution of magnetization as a function 
of the applied electric field. These experiments demonstrate that the strain-mediated 
electric field control of magnetization indeed is a viable technique in technologically 
relevant ferromagnetic thin films at room temperature. 

2. Sample preparation and experimental techniques 

We fabricated ferromagnetic thin film/piezoelectric actuator structures by depositing 
nickel (Ni) thin films onto piezoelectric Pb(Zr x Tii_ x )03-based actuators [24]. Nickel 
was chosen as the ferromagnetic constituent as it is a prototype 3d itinerant ferromagnet 
with a Curie temperature T c = 627 K well above room temperature [25], a high bulk 
saturation magnetization M s = 411 kA/m [26] and sizeable volume magnetostriction [271 
[28] A = |Aioo+|Ain = — 32.9xl0~ 6 with Ai o and A m being the single crystal saturation 
magnetostriction for a magnetic field applied along a crystalline < 100 > or < 111 > 
axis, respectively. The actuators exhibit a hysteretic mechanical stroke of up to 1.3 xlO -3 
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Figure 1. (a) Schematic illustration of the ferromagnetic thin film/piezoelectric 
actuator hybrid, (b), (d) The application of a voltage V p ^ V to the actuator results 
in a deformation of the actuator and the affixed ferromagnetic film. The relaxed 
actuator at V p = V is shown by dotted contours, (c), (e) Schematic free energy 
contours in the film plane. The magnetic easy axis (e.a.) shown by the thick dashed 
line is oriented parallel to the compressive strain and can thus be rotated by 90° by 
changing the polarity of V p . 



along their dominant elongation axis [cf. Fig. [11(a)] if voltages —30 V < V p < +150 V 
are applied. Prior to the deposition of the Ni film, the actuators were mechanically 
polished to a size ofxxt/xz = 3x 2.6 x 2 mm 3 [cf. Fig. CE](a)] to accommodate the 
size of the sample to the restrictions imposed by the ferromagnetic resonance setup. We 
then used electron beam evaporation at a base pressure of 2.0 x 10 _8 mbar to deposit a 
70 nm thick Ni film onto an area of 5 mm 2 on the x-y face of the actuators. To prevent 
oxidation of the Ni film, a 10 nm thick Au film was deposited in-situ on top of the Ni 
layer. The multifunctional hybrid obtained after Ni deposition is sketched schematically 
in Fig. [T](a) together with the definition of the angles that describe the orientation of 
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the magnetization M = (M, G, <3>) and the external magnetic field H = (H, 6, 0) in the 
sample-affixed coordinate system. 

To determine the static magnetic response of the ferromagnetic thin 
film/piezoelectric actuator hybrid we employ superconducting quantum interference de- 
vice (SQUID) magnetometry. The Quantum Design MPMS-XL-7 SQUID magnetome- 
ter is sensitive to the projection m = mH of the total magnetic moment m onto the 
unit vector H = 2. We corrected m for the paramagnetic contribution of the actuator 
and used the Ni film volume V = 3.5 x 10~ 13 m 3 to calculate the respective projection 
M = m/V of the magnetization onto H. All magnetometry data shown in the following 
were recorded at a temperature T = 300 K. 

The magnetic anisotropy of the ferromagnetic thin film/piezoelectric actuator 
hybrid was measured by ferromagnetic resonance (FMR) at room temperature. We 
use a Bruker ESP 300 spin resonance spectrometer with a TE102 cavity operating at a 
constant microwave frequency %w = 9.3 GHz. The sample can be rotated in the FMR 
setup with respect to the external dc magnetic field, so that either 9 or <p [cf. Fig. [D(a)] 
can be adjusted at will. To allow for lock-in detection we use magnetic field modulation 
at a modulation frequency of 100 kHz with a modulation amplitude of /io-^mod = 3.2 mT. 

3. Phenomenology of strain-induced magnetic anisotropy 

The piezoelectric actuator deforms upon the application of a voltage V p 7^ V. Due 
to its elasticity, an elongation (contraction) along one cartesian direction is always 
accompanied by a contraction (elongation) in the two orthogonal directions. Therefore, 
for V p > V, the actuator expands along its dominant elongation axis y and contracts 
along the two orthogonal directions x and z. The Ni film affixed to the x-y face of the 
actuator is hence strained tensilely along y and compressively along x for V p > V [cf. 
Fig. CD(b)]. For V p < V the actuator contracts along y and thus expands along x and 
z and the Ni film thus exhibits a compressive strain along y and a tensile strain along 



To describe the impact of this lattice strain on the Ni magnetization orientation 
we use a magnetic free energy density approach. The free energy F tot is a measure for 
the angular dependence of the magnetic hardness, with maxima in F tot corresponding 
to magnetically hard directions and minima to magnetically easy directions. In 
equilibrium, the magnetization always resides in a local minimum of F tot . Contrary 
to single-crystalline films, the evaporated Ni films are polycrystalline and thus show no 
net crystalline magnetic anisotropy which may compete with strain-induced anisotropies 
and hereby reduce the achievable magnetization orientation effect pjj]. F tot is thus given 



The first term F stat = —fi Q MH (sin 6 sin $ sin 6 sin + cos 9 cos 9 + sin cos $ sin 6 cos (j>) 
in Eq. (pQ) is the Zeemann term and describes the influence of an external magnetic 



x [cf. Fig. Ufa)]. 



by 
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field H on the orientation of M. The uniaxial demagnetization term -F demag = 
^■M 2 sin 2 9 cos 2 $ is the anisotropy caused by the thin-film shape of the sample [29]. 
The last contribution to Eq. (pp) 

Fmagci = (c$ - c™) Msin 2 6 sin 2 $ - 1/3) (2) 

+ £ 2 (cos 2 - 1/3) + £ 3 (sin 2 6 cos 2 $ - 1/3)] 

describes the influence of the lattice strains on the magnetic anisotropy [27]. The strains 
along the x-, y- and z-axis are denoted in Voigt (matrix) notation [30] as £i,£2 and £3, 
respectively. Furthermore, c^ 1 = 2.5 x 10 11 N/m 2 and c^ 1 = 1.6 x 10 n N/m 2 are the 
elastic moduli of Ni [28]. The effects of shear strains (si,i G {4, 5, 6}) average out in our 
polycrystalline film and thus are neglected. 

Using sin 2 sin 2 $ + cos 2 + sin 2 cos 2 $ = 1 and omitting isotropic terms, Eq. ([2j) 
can be rewritten as 



with 



^magel = #u,magcl,y COS 2 + if u , m ag e l,z sin 2 COS 2 $ (3) 

K„i,y=lK^-c^)(s 2 -e 1 ) (4) 
3— ■ ■ 

^u,magel,z = 2^( C 12 — c ll)( e 3 ~~ £ l) ■ 

Due to the elasticity of the actuator and the Ni film, the strains £j are not 
independent of each other. The strains £\ and e 2 in the x-y plane are linked by the 
Poisson ratio v = 0.45 of the actuator [24] according to 

£1 = -v£2 • (5) 

Furthermore, due to the elasticity of the Ni film, the out-of-plane strain £ 3 can be 
expressed as [H] 

c Ni 

£3 = --§(£! + £2) • (6) 
c n 

Assuming a linear expansion of the actuator with lateral dimension L parallel to its 
dominant elongation axis, we can calculate the strain in the Ni film parallel to the 
actuator's dominant elongation axis as 

e = 5L Fp (7) 
2 L 180V ' { ' 

where SL/L = 1.3 x 10~ 3 is the nominal full actuator stroke for the full voltage swing 

-30 V < V p < +150 V. With Eqs. ©, © and © all strains in the Ni film can 

thus directly be calculated for any given voltage V p . This allows to determine the 

equilibrium orientation of the magnetization as a function of V p by minimizing the total 

magnetoelastic free energy density F tot [Eq. (Op)]. 

Due to the negative magnetostriction (A < 0) of Ni and c^ 1 > c^ 2 \ the in-plane easy 

axis of F tot is oriented orthogonal to tensile and parallel to compressive strains in the 

absence of external magnetic fields. Due to the strong uniaxial out-of-plane anisotropy 
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caused by -Fdemag, the in-plane easy axis is the global easy axis. For V p > OV, the 
Ni film exhibits a tensile strain along y (e^ > 0) and a compressive strain along x 
(ei < 0). The easy axis of F tot is thus oriented along the x-direction [cf. Fig. H(c)]. 
Accordingly, for V p < V, the easy axis of F tot is oriented parallel to the compressive 
strain along y and orthogonal to the tensile strain along x [cf. Fig. [0(e)]. We thus 
expect a 90° rotation of the in-plane easy axis of F tot upon changing the polarity of V p . 
Note that in the absence of external magnetic fields there are always two energetically 
equivalent antiparallel orientations of M along the magnetic easy axis. However, for 
simplicity, we only show one of the two resulting possible in-plane orientations of M in 
Figs. [0(b) and [0(d). Nevertheless these panels show that it should be possible to change 
the orientation of M from M| |x to M| |y via the application of appropriate voltages V p 
to the actuator. 

The total magnetic free energy density F tot in Eq. (pp) is experimentally accessible 
by ferromagnetic resonance measurements. The FMR equations of motion [SH [321 EH] 



link the experimentally determined ferromagnetic resonance field fioH ves to the magnetic 
free energy density F tot . In Eq. jEJ), 7 is the gyromagnetic ratio and ui = 27ti/mw with the 
microwave frequency %w The effect of damping is neglected in Eq. ([HI), as damping 
only affects the lineshape which will not be discussed in this article. We note that 
the FMR resonance field HoH ms measured in experiment can be considered as a direct 
indicator of the relative magnetic hardness along the external dc magnetic field direction. 
From Eqs. (DQ) and ([8]) one finds that smaller resonance fields correspond to magnetically 
easier directions and larger resonance fields to magnetically harder directions. 

4. Ferromagnetic resonance 

In this Section, we quantitatively determine the magnetic anisotropy of our ferromag- 
netic thin film/piezoelectric actuator hybrid structure using FMR measurements. Fig- 
ure [2](a) shows four FMR spectra recorded at room temperature with H||y and a volt- 
age V p G {— 30 V, V, 30 V, 90 V} applied to the actuator, respectively. Each spectrum 
shows one strong ferromagnetic resonance which - due to the magnetic field modula- 
tion and lock-in detection - has a lineshape corresponding to the first derivative of a 
Lorentzian line [34] • The resonance field [ioH res of a given FMR line is determined as 
the arithmetic mean of its maximum and minimum and depicted by the full squares 
in Fig. E](a). The figure shows that, with the external magnetic field H||y, HoH res is 
shifted to higher magnetic fields for increasing V Pl while the lineshape is not significantly 
changed. According to Eqs. (pQ) and (jSJ), this implies that the y direction becomes in- 
creasingly harder as V p is increased. 

To determine the evolution of fioH res with V p in more detail, we recorded FMR 
spectra similar to those shown in Fig. [2](a) for V p increasing from —30 V to +90 V 
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Figure 2. (a) FMR spectra recorded with H||y at different voltages V p . An 
increasing resonance field (XoH Ies (solid squares) is observed for increasing V p while 
the lineshape is not significantly altered, (b) The dependence of iioH res on V p is 
qualitatively different for H||y and H||x. Full symbols correspond to increasing V p 
and open symbols to decreasing V p . The solid lines represent the resonance fields 
calculated from magnetoelastic theory (cf. Sec. yielding very good agreement with 
the measurement. For H||y, increasing V p increases ^oH ICS and thus the magnetic 
hardness of this direction, while for H||x, increasing V p decreases fJ,aH Tes and thus the 
magnetic hardness of this direction. The hysteresis of fioH rcs is due to the hysteretic 
stress-strain curve of the actuator. 

(upsweep) and decreasing back to —30 V (downsweep) in steps of AVp = 10 V. For 
H||y, we obtain the resonance fields shown by the squares in Fig. [2(b). Here, full 
squares depict the upsweep and open squares the downsweep of V p . As discussed in the 
context of Fig. E](a), Fig. [2(b) shows that the FMR resonance field for H||y increases 
with increasing V p and decreases with decreasing V p . The small hysteresis between 
up- and downsweep is due to the hysteretic expansion of the actuator plj. Carrying 
out the same series of FMR measurements with H||x yields the resonance fields shown 
by the triangles in Fig. [2](b). For this magnetic field orientation, fioH ves decreases for 
increasing V p and vice versa. In terms of magnetic anisotropy we thus can conclude that 
the x direction becomes easier for increasing V p while at the same time the y direction 
simultaneously becomes harder. 

For a more quantitative discussion we have also plotted the behavior expected from 
Eqs. (jTJ) and (jHJ) as solid lines in Fig.[2](b). These lines show the resonance fields obtained 
by assuming a linear, non-hysteretic voltage-strain relation [cf. Eq. and solving 
Eq. (jHJ) with F to t from Eq. (pp). We use a saturation magnetization M s = 370kA/m 
as determined by SQUID measurements and a g-factor of 2.165 [35]. Considering the 
fact that the actuator expansion saturates at high voltages and the data in Fig. [2](b) 
only show a minor loop of the full actuator swing, the simulation is in full agreement 
to the experimental results. This demonstrates that for Vp < V the y-direction is 
magnetically easier than the x-direction, while for V p > V the x-direction is easier 
than the y-direction. Moreover, at V p = V, both measurement and simulation yield 
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(a) e=o° (b) 0=0° 




Figure 3. (a) The symbols show the FMR resonance field ^oH ros (9) obtained for 
a constant actuator voltage V p = —30 V (open triangles) and V p — +20 V (solid 
circles) as a function of the orientation 9 of H in the sample plane. A uniaxial 
(180° periodic) anisotropy of ^lqH ycs (0) is observed for both voltages. However, the 
easy axis is rotated by 90° as V p is changed from +20 V to —30 V. The full lines 
show the resonance fields simulated using the anisotropy constants from Eq. ([9]). (b) 
Corresponding experiments with H rotated in the y plane (from within the Ni film 
plane to out-of-plane) . Regardless of V p , a strong uniaxial anisotropy with the hard 
axis perpendicular to the film plane is observed. The inset shows that the resonance 
fields for H in the film plane (</> = 90° and (j> = 270°) are still shifted as a function of 
V p in accordance to the data shown in (a) for 8 = 90° and 9 = 270°. 

resonance fields of fioH res m 150 mT for H||y as well as H||x. Thus, at this voltage, 
both orientations x and y are equally easy in our sample. 

To quantitatively determine the full magnetic anisotropy as a function of Vp, we 
recorded FMR traces at constant V p for several different H orientations. The FMR 
resonance fields thus determined in experiment are shown in Fig. [3] together with 
simulations according to Eqs. (pQ) and jBj). The open triangles in Fig. [3] represent 
the resonance fields obtained for V p = —30 V and the full circles those obtained for 
V p = +20 V. Note that in the polar plots in Fig. [3J the distance of fioH res from 
the coordinate origin is an indicator of the magnetic hardness, with easy directions 
corresponding to small distances and hard directions to large distances. 

If H is rotated in the film plane [cf. Fig. [3](a), = 90°], the obtained fi H Tes (9) 
shows minima at 9 = 0° and 9 = 180° for V p = -30 V and at 9 = 90° and 9 = 270° 
for V p = +20 V, respectively. Thus, a clear 180° periodicity of the resonance fields and 
hence a uniaxial magnetic anisotropy is observed for both V p . As the orientations 9 
corresponding to minima of /^oH ves for one voltage coincide with maxima for the other 
voltage, we conclude that the direction of the easy axis is rotated by 90° if V p is changed 
from Vp = —30 V to Vp = +20 V. This is exactly the behaviour expected according to 
Figs.[H(c) and[D(e). 
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If H is rotated from within the film plane (0 = 90°, 6 = 90°) to perpendicular to 
the film plane (4> = 0°, 6 = 90°), we obtain the resonance fields shown in Fig. [3](b). In 
this case, we observe a strong uniaxial anisotropy with a hard axis perpendicular to the 
film plane regardless of V p , stemming from the ^-independent contribution Fdemag to 
F tot . The inset in Fig. [3](b) shows that for H in the film plane, the resonance fields for 
V p = —30 V and V p = +20 V are shifted by approximately 10 mT in full accordance to 
the data in Fig. [3](a). 

As discussed in Section [3J the V p dependence of the resonance field is given by a 
V p dependence of F ma&e \. We thus evaluate the measurement data shown in Figs. [3](a,b) 
using an iterative fitting procedure of i^ u ,magei,y m Eq. (pp) to fulfill Eq. ((ED- We obtain 

#u,magel,y(-30 V)/M s = -3.4 mT (9) 
#u,magel,y(+20V)/M s = 4.4 mT . 

The resonance fields calculated using these anisotropy fields as well as Eqs. (JSj) and (pQ) 
are depicted by the lines in Fig. [3j The good agreement between simulation and 
experiment shows that F tot given by Eq. (pQ) is sufficient to describe the magnetic 
anisotropy of the Ni film. 

In summary, the FMR experiments conclusively demonstrate that it is possible 
to invert the in-plane magnetic anisotropy of our hybrid structure, i.e. to invert the 
sign of -Ku,magei, y solely by changing V p . As the FMR experiment furthermore allows to 
quantitatively determine all contributions to the free energy, Eq. ([TJ), the orientation of 
the magnetization vector M in our sample can be calculated a priori for arbitrary H 
and V p . However, it is not possible to directly measure the magnetization orientation as 
a function of V p with FMR. To demonstrate that the piezo- voltage control of magnetic 
anisotropy indeed allows for a voltage control of M, we now turn to magnetometry. 

5. Magnetometry 

In this Section, we show that it is possible to not only change the magnetic anisotropy 
but to deliberately irreversibly and/or reversibly rotate M simply by changing V p . 
To this end, we need to employ an experimental technique that is directly sensitive 
to the orientation of M rather than to magnetic anisotropies. Here we used SQUID 
magnetometry to record the projection m of the total magnetic moment m onto the 
direction of the external magnetic field H. 

In a first series of experiments we recorded m as a function of the external magnetic 
field magnitude noH at fixed orientations of H and fixed voltages V v at T = 300 K. 
Figure[4](a) shows M = m/V measured with H||y as a function of the external magnetic 
field strength at constant voltage V p = +30 V (full circles) and V p = —30V (open 
triangles). The Ni film shows a rectangular M(H) loop for V p = —30 V, while for 
V p = +30 V the remanent magnetization is lowered by a factor of approximately three. 
According to, e.g., Morrish [29], the rectangular loop for V p = —30 V indicates a 
magnetically easy axis, while the smooth, s-shaped loop for V p = +30 V indicates a 



Voltage controlled inversion of magnetic anisotropy 



10 




-20 -15 -10 -5 5 10 15 20 
H H (mT) 

Figure 4. (a) M (H )-loops recorded with H| |y show a higher remanent magnetization 
for V p = —30 V than for V p = +30 V, thus the y-axis is magnetically easier for 
V p = —30 V than for V p — +30 V. (b) The x-axis is magnetically easier for V p = +30 V 
than for V p = -30 V. 

magnetically harder axis. Thus Fig. [4](a) shows that the y-direction is magnetically 
easier for V p = —30 V and harder for V p = +30 V - as expected from the FMR 
experiments. Changing the orientation of H to H||x yields the M(if )-loops shown 
in Fig. BJb) . Following the same line of argument we can conclude that the x-direction 
is magnetically easier for V p = +30 V and harder for V p = —30 V. Altogether these 
results show that we observe an in-plane anisotropy, the easy axis of which is parallel 
to y for V p = —30 V and parallel to x for V p = +30 V. These observations are fully 
consistent with the FMR results, and corroborate the simple model shown in Fig. [H 

Before discussing further experimental magnetometry results, we note that the free 
energy minima shown in Fig. H](c,e) as well as the in-plane resonance fields in Fig. [3](a) 
are degenerate by 180°. This degeneracy may lead to demagnetization due to domain 
formation if the polarity of V p is repeatedly inverted as there are always two energetically 
equally favorable but opposite directions of M. Thus, to achieve a reversible control 
of M orientation, the degeneracy of F tot needs to be lifted. This can be achieved if a 
unidirectional anisotropy is superimposed on the uniaxial anisotropy. Regarding F tot 
in Eq. (Op), only F stat possesses the desirable unidirectional anisotropy. Hence, a small 
but finite external magnetic field H^0 can be used to lift the 180° degeneracy. This 
approach works for all H orientations, except for H exactly parallel x or y. For H||x 
(with H > 0), the easy axis parallel to x in Fig. [Tj(c) is replaced by an easier positive 
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Figure 5. SQUID M(V p )-loops show the projection M of the magnetization M 
onto the direction of the applied magnetic field H as a function of V p . The symbols 
represent the experimental data and the lines show the simulation of M resulting from 
a minimization of F to t [Eq. {TJ], with £2 determined using a strain gauge to explicitly 
take into account the actuator hysteresis, (a) For H| |y with /io-ff = —3 mT, M exhibits 
an irreversible rotation from A to B followed by a reversible rotation from B to C. (b), 
(c) For H||y and H||x with [iqH = — 5mT, M exhibits a reversible rotation from A 
to B and back to A. 

x direction and a harder negative x direction, thus the degeneracy is lifted for V p > 0. 
However, for V p < [cf. Fig. [T](e)] a magnetic field H||x is orthogonal to the easy axis 
which thus remains degenerate. The same consideration holds for H||y where we expect 
degenerate free energy minima for V p > and a preferred M orientation for V p < 0. 
However, in experiment it is essentially impossible to orient H exactly along x or y; any 
small misorientation between H and x or y is sufficient to lift the degeneracy. 

We now turn to the experimental measurement of the voltage control of M 
orientation. To show that M can be rotated by varying V p alone, we change V p 
at constant external magnetic bias field H. As SQUID magnetometry is limited to 
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Figure 6. Calculated free energy contours in the film plane at the points of the 
M(V p )-loop depicted by capital letters in Fig. [5] (solid lines: V p = +120 V, dotted lines: 
V p = —30 V). To clarify the lifting of the free energy degeneracy by the magnetic field, 
an angle of 10° between the magnetic field and the x or y axis was assumed. The open 
downward-oriented arrows depict the orientation of H during the field preparation at 
7T and the closed downward-oriented arrows depict the orientation of H during the 
actual measurement. Capital letters indicate the equilibrium M orientation at the 
corresponding positions in Fig. [5j (a) Subsequent to the field preparation at —30 V, 
M resides in a local minimum of F to t (point A) at [IqH = — 3mT and rotates to the 
global minimum of F to t at (point B) as V p is increased to V p — +120 V. Sweeping the 
voltage from +120 V to —30 V now results in M following the global minimum of F tot . 
(b), (c) For [IqH = — 5mT, M follows the global minimum of Ftot as V p is changed. 
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recording the projection of M on the direction of the external magnetic field H, it is 
important to choose appropriate H orientations in the experiments. As evident from 
Figs. CD and 31 H||x and H||y are the most interesting orientations of the external 
magnetic field. In view of the circumstance discussed in the previous paragraph, we 
applied H to within 1° of x or y, respectively. We still refer to these orientations of 
H as H||x and H||y for simplicity, but take into account a misalignment of 1° in the 
calculations. As the experimental results and the corresponding simulations will show, 
this slight misalignment is sufficient to lift the degeneracy in the free energy regardless 
of Vp. 

Recording M as a function of V p for two complete voltage cycles —30 V < V p < 
+120 V with H||y and iiqH = — 3mT yields the data points shown in Fig. [5](a). Since 
H||y, M is the projection of M to the y direction in the experiment. Prior to the first 
voltage sweep starting at point A, the voltage was set to V p = —30 V and the sample 
was magnetized to a single domain state by applying fi H = +7T along y which is 
the easy axis for V p = —30 V. The magnetic field was then swept to fioH = —3 mT 
which is close to but still below the coercive field of the rectangular loops [cf. Fig. 0] 
and the acquisition of the M(V P ) data was started. The fact that M is positive at first 
(starting from point A), while the field is applied along the negative y direction, shows 
that M and H are essentially antiparallel at first. Upon increasing V p in steps of +5 V 
to V p = +120 V (point B), M vanishes, which indicates an orthogonal orientation of M 
and H. Upon reducing V p to its initial value of —30 V, M becomes negative (point C), 
evidencing a parallel orientation of M and H. The V p cycle was then repeated once 
more, with M now reversibly varying between its negative value at point C and zero at 
B. Hence the evolution of M is qualitatively different in the first and in the second V p 
cycle, with an irreversible M rotation in the first cycle and a reversible M rotation in 
the second cycle. This behaviour is expected from the magnetic free energy [cf. Eq. (DQ)] 
as the preparation (point A) yields M in a metastable state (M antiparallel H). 

After a renewed preparation with HqH = +7T at V p = —30 V, we repeated the 
experiment with the external magnetic field H||y but with a slightly larger magnitude 
HoH = —5 mT and again recorded M for a complete voltage cycle. The magnetic field 
magnitude of fioH = —5 mT was chosen as it exceeds the coercive field (cf. Fig. [4]) 
while keeping the influence of the Zeemann term in Eq. (pQ) on the total magnetic 
anisotropy comparable to the influence of the strain-induced anisotropics [cf. Eq. ([9])]. 
The experimental data are shown by the symbols in Fig.[5](b). Here a parallel orientation 
of M and H is already observed at point A (V p = —30 V) and M rotates reversibly by 
approximately 70° towards point B (Vp = +120 V) and back when V p is reduced to 
—30 V again (point A). 

To complete the picture, we repeated the M(V P ) experiment again, but now applied 
H||x. The magnetic preparation with fi H = 7T was now performed at V p = +120 V 
to ensure that the preparation field was again applied along an easy axis. The data 
recorded subsequently at \i§H = —5 mT [cf. Fig. E](c)] show a reversible rotation of M 
by approximately 60° between points A, B (Vp = —30 V) and back to A. 
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To quantitatively simulate the evolution of M with V p depicted in Fig. El we again 
took advantage of the fact that the free energy of the Ni film is known from the FMR 
experiments. To account for the hysteresis of the actuator expansion which is responsible 
for hysteresis in the M(V P ) loops in Fig. [5], we hereby used Si as measured using a Vishay 
general purpose strain gauge in the voltage range of —30 V < V p < +120 V. The e 2 
data thus obtained are in accordance with the actuator data sheet. Using e 2 (y p ) as 
well as the material constants given above we obtained in-plane free energy contours for 
each voltage V p . Figure M exemplarily shows such free energy contours at the voltages 
depicted by the capital letters in Fig. [5l To clearly visualize the effect of H misalignment 
with respect to the sample coordinate system on F tot , the plots in Fig. [6] were calculated 
assuming a misalignment in the in-plane orientation of the external magnetic field 9 
of 10°, so that 9 = 10° for H||y and 9 = 100° for H||x. Figure E clearly shows that 
under these conditions the local minima of F tot are non-degenerate for V p = —30 V and 
V p = +120 V. 

To determine the orientation of M, we traced the minimum of the total free energy 
F tot as a function of V p . This was done by minimizing Eq. (pQ) with respect to the 
in-plane orientation 9 of M, whilst assuming that $ = 90° due to the strong shape 
anisotropy. 

To simulate the experimental M(V^,)-loops [cf. Fig. [5], we assume a more realistic 
misalignment of 1° in 9 and numerically minimize -Ftot(@) as a function of V p . For this, 
we set the initial value of antiparallel to the external magnetic field for fi H = —3 mT 
and parallel to the external magnetic field for HqH = — 5mT. Minimizing -F to t(©) 
determines the M orientation which we project onto the y or x axis to yield M. For 
this, the magnitude of M is chosen to give a good fit to the experimental data at points 
A and is assumed to remain constant independent of V p . 

The resulting simulations of M are shown by the solid lines in Fig.El The simulation 
yields the experimentally observed irreversible M rotation during the first voltage sweep 
in Fig. [5](a) as well as the reversible M rotations in the second voltage sweep in Fig. [5](a) 
and in Figs. E](b,c). The simulated total swing of M is in excellent agreement with 
the experimental results for H||y and in good agreement for H||x. The fact that the 
experimental results exhibit a more rounded shape than the simulation is attributed 
to domain formation during the magnetization reorientation, which is neglected in the 
simulation. 

Taken together, the free energy density of Eq. ([TJ with the contributions 
quantitatively determined from FMR and our simple simulation of M(V P ) yield excellent 
agreement with experiment. In particular, we would like to emphasize that for the 
experiment and the simulation the application of V p leads to a rotation of M and not 
to a decay into domains. This is evident from the fact that in Figs. [5](a) and[5](b), a 
large projection of M onto y is accompanied by a small projection onto x and vice 
versa. Combining FMR and SQUID magnetometry, we thus have unambiguously and 
quantitatively demonstrated that M can be rotated reversibly by about 70° at room 
temperature solely via the application of appropriate voltages V p . 
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6. Summary 

The FMR measurements show that the in-plane anisotropy of Ni/piezoactor hybrids 
can be inverted if the polarity of the voltage applied to the actuator is changed. The 
magnetometry results corroborate this result and furthermore show that it is possible 
to irreversibly or reversibly rotate M solely by changing V p . For an irreversible M(V P ) 
rotation, an appropriate preparation of M using a magnetic field sweep is necessary - 
i.e. M must be aligned in a local free energy minimum. It then can be rotated out of 
this minimum by changing V p and the corresponding M orientation change can amount 
to up to 180°. However, this voltage control of M is irreversible in the sense that M can 
not be brought back into the original orientation by changing V p alone. Rather, a second 
appropriate magnetic field sweep is required to align M into its original orientation. In 
contrast, the reversible M(V P ) reorientation requires a preparation of M only once. In 
this case, M is oriented along a global free energy minimum, and can be rotated within 
up to 70° (90° in the ideal case) at will by applying an appropriate V p . The M(V P )- 
loops simulated by minimizing the single-domain free energy [Eq. (DO)] are in excellent 
agreement with experiment, showing that the M orientation at vanishing magnetic field 
strengths can still be accurately calculated from the free energy determined from FMR. 
Finally, we would like to point out that the hysteretic expansion/contraction of the 
actuator visible as a hysteresis in M{V P ) [cf. Fig. [5] in particular also leads to distinctly 
different M{V P = 0), depending on the V p history. Thus, our data also demonstrate that 
a remanent M control is possible in our ferromagnetic thin film/piezoelectric actuator 
hybrids. 
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